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We analyze in detail all the possible fixed points of the effective Hamiltonian of a non-magnetic 
impurity hopping between two sites in a metal obtained by Moustakas and Fisher(MF). We find a 
line of non-fermi liquid fixed points which continuously interpolates between the 2-channel Kondo 
fixed point(2CK) and the one channel, two impurity Kondo (21K) fixed point. There is one relevant 
direction with scaling dimension 1/2 and one leading irrelevant operator with dimension 3/2. There 
is also one marginal operator in the spin sector moving along this line. The marginal operator, 
combined with the leading irrelevant operator, will generate the relevant operator. For the general 
position on this line, the leading low temperature exponents of the specific heat, the hopping sus- 
ceptibility and the electron conductivity Ci mp , Ximpy °~(T) are the same as those of the 2CK, but the 
finite size spectrum depend on the position on the line. No universal ratios can be formed from the 
amplitudes of the three quantities except at the 2CK point on this line where the universal ratios can 
be formed. At the 2IK point on this line, <r(T) ~ 2a u (l + aT 3 ' 2 ), no universal ratio can be formed 
either. The additional non-fermi liquid fixed point found by MF has the same symmetry as the 21K, 
it has two relevant directions with scaling dimension 1/2, therefore also unstable. The leading low 
temperature behaviors are Cimp ~ T, Ximp ~ log T, o~(T) ~ 2er u (l + aT 3 / 2 ), no universal ratios can 
be formed. The system is shown to flow to a line of fermi-liquid fixed points which continuously 
interpolates between the non-interacting fixed point and the 2 channel spin-flavor Kondo fixed point 
(2CSFK) discussed by the author previously. The effect of particle-hole symmetry breaking is dis- 
cussed. The effective Hamiltonian in the external magnetic field is analysed. The scaling functions 
for the physical measurable quantities are derived in the different regimes; their predictions for 
the experiments are given. Finally the implications are given for a non-magnetic impurity hopping 
around three sites with triangular symmetry discussed by MF. 



I. INTRODUCTION 



The experimental realization of overscreened multichannel Kondo model has_been vigorously searched since the 
discovery of its non-fermi liquid behavior (NFL) by Nozieres and Blandin (NB)El. D. L. Coxa pointed out that the 
NFL behaviors in heavy fermion systems like Y\^ X \J x Pd^ may be explained by the 2 channel quadrupolar Kondo 
effects, but the observed electrical conductivity of suck systems is linear in T in contrast to VT behavior of the 
2-channel Kondo model (2CK). Vladar and Zawadowskid suggested that a non-magnetic impurity tunneling between 
two sites in a metal can be mapped to the 2CK in which the roles of channels and spins in the original formulation 
are interchanged. Ralph et alu proposed that the conductance signals observed in ballistic metal point contacts may 
be due to the 2-channel Kondo scattering from 2-level tunneling systems, the conductance exponent 1/2 and the 
magnetic field dependence observed in such device is iiftked in consistent with that predicted by Affleck-Ludwigis 
(AL) Conformal Field Theory (CFT) solution of the 2CKBEI however the alternative interpretation was also proposed!! 

Moustakas and Fishertj reexamined the problem of the electron assisted tunneling of a heavy particle between two 
sites in a metal. In addition to barg hopping ( Ao term in Eq. |l| ) and one electron assisted hopping term ( Ai term 
in Eq. ^ ) found previously in Ren3, they found that an extra two electrons assisted hopping term ( A2 term in Eq. 
[l] ) also plays an important role. Treating all these important processes carefully, Ahey concluded that more than 
four channels (including spin) are needed in order to localize the impurity. In Ref.tj, they wrote down an effective 
Hamiltonian which includes all the important processes and employed Emery-Kivelson (EK)'s Abelian Bosonization 
solution of the 2CKEj to investigate the full phase diagram of this Hamiltonian, they found the two electron assisted 
hopping term plays a similar role to the bare hopping term. However, they overlooked the important fact thai 
the canonical transformation operator U — e lS * s in EK's solution is a boundary condition changing operatoiO, 
therefore their analysis of the symmetry of the fixed points and the operator contents near these fixed points are 
not complete. They didn't calculate the electron conductivity which is the most important experimental measurable 
quantity. Furthermore the nature of the stable Fermi liquid fixed point was also left unexploited. 
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Affleck and Ludwig (AL) lilHlSEHl, using Conformal Field Theory, pointed out that for any general quantum impurity 
problem, the impurity degree of freedoms completely disappear from the description of the low temperature fixed point 
and leave behind conformally invariant boundary conditions. CFT can also be used to classify all the possible boundary 
operators near any low temperature fixed points and calculate any correlation functions. For 4 pieces of bulk fermions 
which correspondlo 8 pieces of Majorana fermions, the non-interacting theory possesses 50(8) symmetry, Maldacena 
and Ludwig (ML)t3 showed that finding the symmetry of the fixed points is exactly equivalent to finding the boundary 
conditions of the 8 Majorana fermions at the fixed points, the boundary conditions turned out to be linear in the 
basis which separates charge, spin and flavor. ML reduced the descriptions of the fixed points as free chiral bosons 
plus different linear boundary conditions. The linear boundary conditionSpCan also be transformed into the boundary 
conditions in the original fermion basis by the triality transformation Eq.g|!£l. The boundary conditions in the original 
fermion basis only fall into two classes: NFL fixed points where the original fermions are scattered into spinors; fermi 
liquid(F~L) fixed points where the original fermions only suffer just phase shifts at the boundary. 

The important step in the CFT approach developed by AL is the identification of the fusion rules at various fixed 
points. Although the fusion rule is simple in the multichannel Kondo model, it is usually very difficult to identify in 
more complicated models like the one discussed in this paper. 

Recently, using EK's Abelian Bosonization approach to the 2CK, the author developed a simple and powerful 
method to study certain class of quantum impurity models with 4 pieces of bulk fermions. the method can identify 
very quickly all the possible boundary fixed points and their maximum symmetry, therefore circumvent the difficult 
tasks to identify the fusion rules at, different fixed point or line of fixed points, it can also demonstrate the physical 
picture at the boundary explicitly!! 3 ! _ 

In this paper, using the method developed in Ref.t 3 ] and paying the special attention to boundary condition changing 
nature of U = e lS * s , we investigate the full phase diagram of the present problem pagain. In Sec. II, we Abelian 
bosonize the effective Hamiltonian. By using the Operator Product Expansion (OPE)liil, we get the Renormalization 
Group(RG) flow equations near the weak coupling line of fixed points , therefore identify the two independent crossover 
scales. In the following sections, we analyze all the possible fixed points of the bosonized Hamiltonian. In Sec. Ill, 
we find a line of NFL fixed points which-,continuously interpolates between the 2CK fixed point and the one channel 
two impurity Kondo (2IK) fixed pointed, its symmetry is U(l) x 0(1) x 0(5)113. This line of NFL fixed points is 
unstable, it has one relevant direction with scaling dimension 1/2. It also has one marginal operator in the spin sector 
which is responsible for this line. The OPE of this marginal operator and the leading irrelevant operator will always 
generate the relevant operator. For the general position on the line, although the leading exponents of the specific 
heat, hopping susceptibility and the electron conductivity Oi mp , x^ mp , er(T) are the same as those of the 2CK, the 
finite size spectrum depend on the position on the line; no universal relations can be found among the amplitudes 
of the three quantities. Only at the 2CK point on the line, universal ratios can be formed. However, at the 2IK 
point, the coefficient of y/T vanishes; we find two dimension 5/2 operators which lead to a(T) ~ 2er u (l + T 3 / 2 ), 
no universal ratio can be formed either. In Sec. IV, the additional NFL fixed point found by MS is shown to have 
the symmetry 0(7) x 0(1), therefore is the same fixed point as the 2IK. This fixed point is also unstable, it has 
two relevant directions with scaling dimension 1/2. Because the leading irrelevant operators near this NFL fixed 
point are first order Virasoro descendant with scaling dimension 3/2 which can be written as a total imaginary time 
derivative, therefore, can be droped. The subleading irrelevant operators with dimension 2 give Ci mp ~ T. However, 
because the 'orbital field' in Eq couples to a non-conserved current, Ximp ~ l°g^- We also find two dimension 
5/2 irrelevant operators, one of them contributes to the leading low temperature conductivity a(T) ~ 2cr u (l + T 3 / 2 ). 
No universal ratios can be formed near this NFL fixed point. In Sec. V, we find the system flows to a stable line of 
FL fixed points which continuously interpolates between the i»pn- interacting fixed point and the 2 channel spin-flavor 
Kondo (2CSFK) fixed point discussed by the author in RefrJ, its symmetry is U(l) x 0(6) ~ U(A). Along this line 
of fixed points , the electron fields of the even and the odd parity under interchanging the two sites suffer opposite 
continuously- changing phase shifts. We also discuss the effect of the marginal operator in the charge sector due to 
the P-H symmetry breaking and compare it with the magical operator in the spin sector. In Sec. VI, we analyse the 
effective Hamiltonian in the external magnetic field which break the channel symmetry. In Sec. VII, all the scaling 
functions for the physical measurable quantities including the real spin susceptibility are derived in different regimes. 
In the final section, the relevance of the results of this paper to the experiments are examined, some implications on 
the non-magnetic impurity hopping around three sites with triangular symmetry are also given. In Appendix A, the 
finite size spectrum of one complex fermion is listed. In Appendix B, the boundary conditions in the original fermion 
basis are derived by both Bosonization method and 7 matrix method. In Appendix C, the results on the additional 
NFL fixed point in Sec. IV. are rederived in a different basis. 
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II. BOSONIZATION OF THE EFFECTIVE HAMILTONIAN AND THE WEAK COUPLING ANALYSIS 



We start from the following effective Hamiltonian for a non-magnetic impurity hopping between two sites in a metal 
first obtained by MFOJ3 : 

H = H + V 1 (4 (T ^ 1(7 +4^2a) 

+ 4h i^r + -y zZ + ^TTT^l^Mi^n) + h - c - + ■■■ (!) 

c a 

Here the two sites 1,2 ( the two real spin directions tj-L ) play the r °l e of the two spin directions f, j ( the two 
channels 1,2 ) in the magnetic Kondo model. All the couplings have been made to be dimensionless. As emphasized 
by ME- even initially A 1? A 2 maybe negligible, they will be generated at lower energy scales, • • • stands for irrelevant 
termal3. In the following, we use the notation of the magnetic Kondo model and rewrite the above Hamiltonian as: 

H = H + ViJ c (0) + 2V 2 J X (0) + W 3 S z J z (0) + A 1 (J X (Q)S X + J y {0)S y ) 

+ —S x + A 2 2 7 rr c (5-^ 1 t T V'uV'2T^2| + h.c.) + h( [ dxJ z (x) + S z ) (2) 

TTT-c 11 J 

where S + = d\d 2 ,S~ = d\di,S z — k(d\di — d\d 2 ). We also add a 'uniform field' which corresponds to strain or 
pressure in the real experiment of Ref.Q. The real magnetic field will break the channel symmetry, its effect will be 
discussed in Sec. VI. 

Hamiltonian Eq. || has the global Z 2 x SU/(2) x U c (l) symmetry and Time reversal symmetry. Under the,-^ 2 
symmetry in the spin sector : ip^ < — > -0^,S ,a: — > S X ,S V — > —S V ,S Z — > — S z . Under Time reversal symmetrylia : 
i — > —i,ipL — * i/jr,S x — > S X ,S V — > — S V ,S Z — ► 5 Z . The potential scattering term Vi is the only term which breaks 
P-H symmetry: ifii a (x) — > e a /3ipj /3 (x). In contrast to the 2CK, we only have Z 2 symmetry in the spin sector, the total 
spin current to which h couples is not conserved r-. 

In the following, we closely follow the notations of Emery-Kivelsontd. Abelian-bosonizing the four bulk Dirac 
fermions separately : 

iMz) = J^ e -^(z) (3) 
V2ttt c 

Where &i a (x) are the real chiral bosons satisfying the commutation relations 

[$ia(x),®jp(y)] = Si 3 S a f3iTTsgn(x - y) (4) 

The cocycle factors have been chosen as: P 1T = P n = e" Nl \P 2] = P24 = e lw( - Nl ^ +N ^ +N2 i\ 
It is convenient to introduce the following charge, spin, flavor, spin-flavor bosons: 




*1T + *U + $2T + $21 

1 I $ 1T - $ u + $ 2T - $ 2i 

2 [ $ 1T + $ u - $ 2T - $ 2i 



(5) 



Following the four standard steps in EK solution: step 1 : writing the Hamiltonian in terms of the chiral bosons 
Eq.[|. step 2: making the canonical transformation U — exp(— iVi$ c (0) + iS z Q s (0)). step 3 shift the spin boson by 
d x ^s — * dx^s + jjjr"- s t e P 4 : making the following refermionization: 

S x = -^=e inN 'f , S y = -=e™ N °s , S z = ~iab = d*d - 1 
V2 ' \/2 2 

V2 V^TT^c 

^ = -4(a s ,i ~ = -^e^W'te- 1 *' (6) 
V2 V 27rT c 

Note %l> s ,i(x) defined above contains the impurity operator e l7Td ^ d i n order to satisfy the anti-commutation relations 
with the other fermions. 
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The transformed Hamiltonian W — U HU 1 can be written in terms of the Majorana fermions as: 

H' = H + 2yaba s . l (0)b sf (0) + 2qaba s , l {0)b s . l {0) - iabq— 

- l -^ab s} {Q) + i-£±=ba Bti (0) + i-^=ab s>i (0) (7) 
where y = 2V 2 ,q = ±(V 3 - A± = A ± A 2 . 

As observed by MF, the above equation clearly indicate that the two electron assisted hopping term plays a similar 
role to the bare hopping term. From the OPEci of the various operators in Eq. [j], the R. G. flow equations near the 
weak coupling fixed point q = isc3 



dA+ 1 A 

^r = 2 A + + 2yAl 

dA_ 1 



dl 2 
dAi 1 



dl 2 
dy 

dl 
dq 

~dl 



A_ 

Ai + 2yA+ 
AiA + 

A+A_ (8) 



MF got the same RG flow equations ( Eq.23 in Ref.B ) near the weak coupling line of fixed points by using Anderson- 
Yuval Coulomb gas picture. Eq. || shows that A + , A_, Ai have the same scaling dimension 1/2 at q = 0, so are 
equally important. 

Define b S} i(x),b s f,i(x) (similarly for q,y) 



Where Ajf = y'Af + Al 
Eq.fj] can be rewritten as: 



b, ti (x) \ _ _•_ ( Ai A_ \ / b.,i(x) 
b sU (x) J A K \ -A- Ai J\ b sf (x) 



H' = H + 2qaba Sli {0)b Sli (0) + 2yaba Sti (0)b sfti (0) 
A K _~ . A+ - h 



(9) 



-i- 



V27TT C \/2lTT c VF 

The R. G. flow equations which are equivalent to Eq.|| are (we set h = 

dA+ 1 



a,b s f,i(0) + i =ba S; j(0) — iabq — (10) 



dl 2 
dA K 1 



dl 2 
d6 _ A + 

dl ~ q A^ 
dy 

dl 
dq 
~dl 



A+ + 2yA x 
A K + 2yA^ 



A K A+ 

(11) 



Where the angle 9 is defined by cos 9 — A A2 Al , sin 6 = 2A . 2 Al . 

"if "K 

The crossover scale from the weak coupling fixed point q — to a given point on the line of NFL fixed points to be 
discussed in section III is given by Tk\ ~ D(Ak) 2 , to the additional NFL fixed point to be discussed in section IV. 
is given by T K2 ~ fl4A+) 2 - 

As emphasized inEJO, the canonical transformation U is a boundary condition changing operator, the transformed 
field tp' s (x) is related to the original field -0s (x) byE-J: 
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ij)' a {x) = U'^s^U = e lwdU e lTrS * S9nx ij s (x) = -isgnx^ s {x) (12) 

As expected, the impurity spin S z drops out in the prefactor of the above equation. 
The above Eq. can be written out explicitly in terms of the Majorana fermions 

af(0) = -6f(0), ^(0)=of(0) 

af (0) = 6f (0), bf (0) = -of (0) (13) 
We hnd the physical picture can be more easily demonstrated in the corresponding action: 

S = S + -£ [ dra(r)^- - f dra(T)b sU (0,r) 

Z J OT y/2lTT c J 

+ |/^)^ + ^/^K i (0,r) 

+ 2q [ dTa(T)b(T)a s ,i(0,T)b s ,i(0,T) + 2y [ dTo(T)6(T)o.,i(0,T)6 a/li (0,r) (14) 



When performing the RG analysis of the action S, we keepc3 1: 72 = 1, Ax fixed, 2: 71 = 1, A + fixed, 3: Ak, A + 
fixed. 

We will identify all the possible fixed points or the line of fixed points and derive the R. G. flow equations near 
these fixed points in the following sections respectively. 

III. THE LINE OF NFL FIXED POINTS 

If q = y = A + = 0, this fixed point is located at 71 = 0,72 = 1 where b decouples, but a loses its kinetic energy 
and becomes a Grassmann Lagrangian multiplier, integrating a out leads to the boundary conditions : 

b L sf ,M = -bf u (0) (15) 

We also have the trivial boundary conditions 

< i (0) = < i (0), 6^(0) = &«(0) (16) 

By using Eqs.p|JT^, we find the boundary conditions in i/'El 

•?<•>-«?«». <"» 

By using Eq.[l3], we get the corresponding boundary conditions in H 

b L s (0) = -6f (0), tf(0) = e i2 V« (0), = 25 (18) 

where Va(0) = a s (0) - ib sf (0). 

It is evident that at the fixed point, the impurity degree of freedoms totally disappear and leave behind the 
conformally invariant boundar y cpr j di jl| iQ ns Eq.|l8. These are NFL boundary conditions. In contrast to the boundary 
conditions discussed previously E3llj'HE3, although the boundary conditions are still linear in the-basis which separates 
charge, spin and flavor, they are not in any four of the Cartan subalgebras of 5*0(8) groupEj, therefore cannot be 
expressed in terms of the chiral bosons in Eq.|5|. 

However, Eq.|l8| indicates that it is more convenient to regroup the Majorana fermions a sMM 

C = a, - ib sf = e~ 1 *" 

C/ = a af + ib s = e-**?/ (19) 
The boundary condition Eq.[l8| can be expressed in terms of the new bosons 

Kl = K« + °i ^f.L = ~*sf,R ( 20 ) 
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As pointed out in RefE in the basis of EqJ^, the physical fermion fields transform as the spinor representation 
of SO(8), therefore in order to find the corresponding boundary conditions in the physical fermion basis, we have to 
find the 16 x 16 dimensional spinor representation of boundary conditions Eq. [jj. The derivation of the boundary 
conditions is given in Appendix B, the results are found to be : 



j,L _ ±18/2 oR 



„±i0 



(21) 



Where the new fields are defined by 



Ipi: 



V2 



Si± = -j= (s^ ± s n) 



(22) 



It can be checked explicitly the above boundary conditions satisfy all the symmetry requirement ( namely Z2 x 
SUf(2) x U c (l), Time Reversal and P-H symmetry). Fermions with the even and the odd parity under Z2 symmetry 
are scattered into the collective excitations of the corresponding parity which fit into the S spinor representation of 
SO(8), therefore the one particle S matrix and the residual electrical conductivity are the same with those of the 2CK. 
This is a line of NFL fixed points with g = v2 and the symmetry O(l) x U(l) x 0(5) which interpolates continuously 
between the 2CK fixed point and the 2IK fixed point. If 9 = tt, namely A_ = 0, the fixed point symmetry is enlarged 
to 0(3) x 0(5) which is the fixed point symmetry of the 2CK. If 6 = 0, namely Ai = 0, the fixed point symmetry is 
enlarged to 0(1) x 0(7) which is the fixed point symmetry of the 2IK(b) (Fig.|j). The finite size —l<x<l spectrum 
( in terms of unit : ^ E - ) can be obtained from the free fermion spectrum with both periodic (R) and anti-periodic (NS) 
boundary conditions by twisting the three Majorana fermions in the spin sector: changing the boundary condition of 
the Majorana fermion b s from NS sector to R sector or vice versa, twisting the complex fermion ip s by a continuous 
angle 9 = 26. The finite size spectrum of one complex fermion is derived in Appendix A. The complete finite size 
spectrum of this NFL fix line is listed in Table | if < S < | and in Table |ll|if--| < 6 < 0. The ground state 

energy is E 
Ei — E a 



1 + & 

2 ' 



jg + 5(7 ) 2 with degeneracy d — 2, the first excited energy is E\ — Eq 



Mi (if-f <<5< f ) or 



(if -f < 5 < -f) with d = 4. 



If 5 = ±j(5 = 0), the finite size spectrum of the 2CK ( the 2IK) is recovered. The finite size spectrum of the 2CK 



and the 2IK are listed in Tables |III| and |TV| respectively. _ . 
The local correlation functions at this line of NFL fixed points aretj 



(a(r)a(0)) = -, 

T 



(6^(0,^6^(0,0)) = ^ 



(23) 



We can also read the scaling dimension of the various fields [6] = 0, [a] = 1/2, [a Sy i] — [b s ,i] — 1/2, [6;,^] = 3/2. 



2IK(b) ( 0(l)xO(7) ) 
0=0 



U(l)xO(l)xO(5) 
NFL Fixed line 



2CK ( 0(3)xO(5) ) 

e=jt 




9=0 
Free ( 0(8) ) 



2CSFK ( 0(2)xO(6) ) 



2IK(a) ( 0(l)xO(7) ) 
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FIG. 1. Phase diagram of a non magnetic impurity hopping between two sites in a metal. The line of NFL fixed points 
has the symmetry O(l) x 17(1) x 0(5) with g = \/2 which interpolates continuously between the 2CK fixed point and the 2IK 
fixed point. If 6 — n, the fixed point symmetry is enlarged to 0(3) x 0(5) which is the fixed point symmetry of the 2CK. 
If 6 — 0, the fixed point symmetry is enlarged to O(l) x 0(7) which is the fixed point symmetry of the 2IK(b). This line of 
NFL fixed points is unstable, there is one relevant operator with dimension 1/2 which drives the system to the line of FL fixed 
points . There is also a marginal operator along the line. This line of FL fixed points has the symmetry (7(1) x 0(6) ~ (7(4) 
with g = 1 which interpolates continuously between the non-interacting fixed point and the 2CSFK fixed point. If 8 — 0, the 
fixed point symmetry is enlarged to 0(8) which is the fixed point symmetry of the non-interacting electrons. If 6 = n, the 
fixed point symmetry is enlarged to 0(2) x 0(6) which is the the fixed point symmetry of the 2CSFK. This line of FL fixed 
points is stable. There is a marginal operator along this line. The additional NFL fixed point with g — v2 has the symmetry 
O(l) x 0(7) which is the fixed point symmetry of the 2IK(a). This additional NFL fixed point is also unstable. There are two 
relevant terms with scaling dimension 1/2, any linear combination of the two terms will drive the system to a given point of 
the line of FL fixed points . See the text for the detailed discussions on the physical properties of these fixed points or line of 
fixed points s and the crossovers between them. 

As shown in Ref.lll, at the line of fixed points , the impurity degree of freedoms completely disappear: b decouples 
and a turns into the non-interacting scaling field at the fixed point 

S(r)~6 5/>i (0,T) (24) 

The corresponding two scaling fields in H is 

-^(-A_a s + Ajbsf) (25) 
Following RefS, we find the impurity spin turns into 

S x (t) ~ ~i(bb s + ^^a s b sf ) H 



A 



S v {t) ~ i(ba s + — !— (-A_a s + Ai& s /)6 S ) H 

S z (T)~ib^-(-A_a s + A 1 b sf ) + --- (26) 



A" 



Where • • • stands for higher dimension operatorscS and = y 1 °° s9 , = y 1+c ° s . 
The impurity spin-spin correlation function (S a (T)S a (0)) ~ 1/r. 

From Eq.0, it is easy to see that even y term itself is irrelevant near this line of fixed points , but Eq. || shows that 
it, when combined with Ao term, will generate Ai, A2 terms which must be taken into account at this line of fixed 
points . y term is " dangerous" irrelevant, similar " dangerous" irrelevant term occurred in the two channel flavor 
anisotropic Kondo modeES. This line of NFL fixed points-is unstable. A + term in Eq. |l4] has scaling dimension 1/2, 
therefore is relevant, this term was first discovered by MFH. The OPE of a s ^ with itself will generate the dimension 2 
energy momentum tensor of this Majorana fermion a S i(0, t) 9as ^^ T ) ■ the OPE of this energy momentum tensor with 

a St i will generate a first order descendant field of this primary field with dimension 3/2 i_ia Si i(0, r) = 8as g^°' T ' > izl. q 
term is the leading irrelevant operator with scaling dimension 3/2, therefore contributes to 

C lmp ~ {qfTlogT (27) 



Where q = J^=f^-q + x P^f^y 



It is important to see there is a marginal operator (9$"(0) in the spin space which changes the angle 9 in Eq,20 
This operator is very different from the marginal operator 9<f' c (0) in the charge sector which changes the angle 9 p h in 
Eq.^5|. Combined with the leading irrelevant operator, it will always generate the dimension 1/2 relevant operator. 
This indicates that the existence of the line of NFL fixed points and the existence of one relevant operator are 
intimately connected. 

However, from Eqsjl^, |2^, only q term contributes to the impurity susceptibility 

xL P ~VlogT (28) 

As shown in RefS, the Wilson Ration R = 8/3 is universal for any general spin anisotropic 2CK. However, near 
this line of NFL fixed points , R is not universal. 
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71 term has dimension 2. y term has scaling dimension 5/2, in H 1 , it can be written as 

y : a(r)^l : a s>i (0, r) ~: 6,^(0, t) ^ ( °' t) : a.,i(0, r) (29) 



1 — cos ( 



Where y = -J±±f^.q- 

This operator content is totally consistent with the following CFT analysis on the special point = ir: if 8 = ir, the 
fixed point symmetry is Q/3) x 0(5), the CFT analysis of Ref.Ea can be borrowed. Under the restriction of Z 2 and 
Time- Reversal symmetrytLa, it is easy to see there is only one relevant operator 1 with scaling dimension 1/2 which 
can be identified as a s ^, two leading irrelevant operators with scaling dimension 3/2, one is the spin operator T° 

which is Virasoro primary and can be identified as q term, another is the spin 1 operator L-icj) 1 = which is the 
first order Virasoro descendant. I— . 
By Eq| and Eq.|f, the spin-0 and spin-1 leading irrelevant operators in H' becomeEa 

qa(T)b{T)a s ,i(0,T)b Sii {0,T) ~ & s/>i (0 ) r)a s , i (0,r)& s , i (0 ) r) 
= & Sl< (0, T)b sf (0, r)a s>i (0, r) = 6' S (Q, r)6 s/ (0, r)a' s (0, r) 
9a s ,i(0,r) cV s (0,r) 



<9r dr 



(30) 



By Eq.13, the corresponding spin-0 and spin-1 operators in H are 



qb s (0,r)(a s (0,r)b sf (0,r)), ^ r) (31) 

Obviously, only the coefficient g depend on the angle which specifies the position on the line of NFL fixed points . 
They can be written in terms of the new bosons introduced in Eq.|l9|, 



cos<F s V(0,T)d<^(0,T), jLcosd^r) (32) 



It is easy to see,that the marginal operator <9<I>"(0) makes no contribution to the Green function. 
Following Rcf.fL2l, the first order correction to the single particle L-R Green functions due to the spin-0 operator can 
be calculated ( x\ > 0, x 2 < ) 

(ip 1+ {x 1 ,Ti)ipl + (x2,T 2 )) ~q J dr{ij 1+ {x l ,T 1 )cos§ n sf (0,T)d<S> n s (0,T)tlj\ + {x 2 ,T2)) 

X ( e -5*/( a: i^i) e 5 <I '/( :I: 2,r2)^ e -^*" / (xi,ri) e i*™ / (0,r) e -±$^ / (2 :2 ,r 2 )^ 

~qe w / 2 ( Zl -z 2 )- 3 / 2 (33) 

where z\ = t\ + ix\ is in the upper half plane, z 2 = t 2 + ix 2 is in the lower half plane. 

The first order correction to the single particle L-R Green functions due to the spin-1 operator vanishes because the 
spin-1 leading irrelevant operator can be written as a total derivative and the three point functions are the periodic 
function of the imaginary time. I— . 

The bosonized form of the y term in Eq.|2^ in H isE3 

y[cos9 : cos 2^(0, r) : ~ : (c>^(0, r)) 2 : 

-sin0(: sin2$' s l (0, r) : -d 2 <I>' s l (0, t))] cos $^(0, r) (34) 
The first order correction due to this dimension 5/2 operator is 

( ipi+ (a>i,nM+ (a?2, t 2 )) ~ 
y J drie-i^^^i-^ : (9$™(0, r)) 2 : + sintfd 2 $™(0, T )] e i^^^+^) 
x ( e -4*?/(*i.n} e »: / (o 1 r) e -4»;>(» a ,7i)j ^y e *e/2( lCl +C 2 sin0)(zi - z 2 y b ' 2 (35) 
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Where C\, C 2 are real numbers. 

From Eq.|32|, it is easy to identify another dimension 5/2 operator 



(a$«( ,r)) 2 :cos$;V(0,r) (36) 



The contribution from this operator can be similarly calculated. 
By using the following OPE: 



. e -i*(*0 gi*^) := (zi _ Z2 )-i/4 _ * (j8l _ Z2) 3/4 . 9$(Z2) 



- ^) 7/4 : 5 2 $(z 2 ) : - z 2 ) 7 ' A : (5$(z 2 )) 2 : + • ■ • (37) 

4 o 



We get the three point functions 



(z 1 -z)(^ 1 -z 2 )-3/ 4 (z-z 2 ) 



(zi - z)(zi - z 2 ) J / 4 (z - z 2 ) 
By Conformal Ward identity, we can write the three point function with the energy momentum tensor 

< e -i*<*0 TWe i»C«.)) = IZ*^ (39) 

In Ref.i, AL found that in the multi-channel Kondo model, the electron self-energy has both real and imaginary 
parts which are non-analytic function of the frequency to. In the presence of P-H symmetry, the imaginary part is 
even function of lu, the real part is odd function of lu, because the two parts are related by Kramers-Kronig relation. 
Only the part of self-energy which is both imaginary and even function of lu contributes to the electron conductivity. 
The factor i will interchange real and imaginary part. In evaluating Eq.^, we used the important fact that the two 
three point functions in Eqs. |38| , |3"9| d iffer by the factor i. 

By conformal transforming Eq.|33| to finite temperature, we get the leading term of the low temperature conductivity 
from channel one and parity + fcrmions 

oi+(T)~ — (1-gsm-VT), a u = (40) 

I I 6ni 

Where p F is_the density of states per spin per channel at the fermi energy, is the impurity density. 

Similarly^, we get the leading term of the low temperature conductivity from channel one and parity — fermions 

^(^-^(l-qsiJ-VT) (41) 

Because of the global SU{2)j symmetry in the flavor sector, the same equations hold in channel 2. 
Even the symmetry in the spin sector is 0(1) x U(l) instead of 0(3) of the 2CK, the 2 channel and 2 parity fermions 
do make the same leading contribution to the total conductivity 

a(T)~2a u (l-qsm-VT) (42) 



For 9 — 7T, namely at the 2CK, q = q, then two universal ratios can be formed from Eqs. [27[|2^ . |42| . 

For 9 — 0, namely at the 2IK, q = y, the coefficient of y/T vanishes. It is evident that the 2nd order correction ( 
actually any even order correction ) to the Green function vanishes, the 3rd order correction will lead to T 3 / 2 , but 
the coefficients still vanishes due to ~ sin 9/2 = 0, because Mdd order corrections have the same i factor. 

By conformal transforming Eq.^35| to finite temperatureEj, we get the next-leading term of the low temperature 
electrical conductivity 

0-1+ - £(cos Q -C x + sin - sin 6»C 2 )T 3 / 2 (43) 



Putting 9 = ( then y = —q ) in the above equation and adding the contribution from the operator in Eq.36 
get the leading term at the 2IK fixed point: 



we 
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(44) 



It is evident that even at the 2IK point, no universal ratios can be formed. 

For general 9, the leading low temperature behaviors of the three physical measurable quantities are given by 
Eqs.^?],^8],[l2], no universal ratios can be formed. 

The potential scattering term Vi is a marginal operator which causes a phase shift in the charge sector: 

$c,l = $ c ,fl + ph (45) 
The symmetry of the fixed point is reduced to O(l) x U(l) x 0(3) x C/(l), Eqs.[h| [H] become: 

CTl 4T)~^(l-gsin^Vr) (46) 

It is easy to see that in the presence of P-H symmetry breaking, in contrast to the 2CK, the different parity fermions 
do make different contributions to the conductivity, Eq^ becomes 

a (T) ~ 2cr u (l — qsin ^ cos ^-VT) (47) 



Eq.[43] becomes 



~, U + V p h„ . V + V p h . „„ ^3/2 
cti+ ~ y{cos 1 2~ sm ® C 2) T 



The total conductivity becomes 



y( cos LJpEltCx + sin LJp* sin 6C 2 )T 3 ' 2 (48) 



~ j/ cos (cos -Ci + sin - sin 6C 2 )T 3/2 (49) 
The total conductivity at the 2IK becomes 

a(T) ~2<7„(1 + cos ^T 3 / 2 ) (50) 

As shown by AL@, P-H symmetry breaking does not change the leading results of the specific heat and susceptibility. 

In this section, we discussed two marginal operators, one in spin space <9<&™(0), another in the charge space <9$ c (0). 
Both make no contributions to the conductivity and make only subleading contributions to the thermodynamic 
quantities. However, d$™(0) is much more important. Combined with the leading irrelevant operator in the spin 
sector, it will generate the dimension 1/2 relevant operator which will always make the line of NFL fixed points 
unstable. Furthermore, as shown by Eqs.^,^8], even the coefficients of the leading terms of the thermodynamic 
quantities depend on the position on the line caused by <9$"(0). 



IV. ADDITIONAL NFL FIXED POINT 



If q = y = Ak = 0, this fixed point is located at 71 = 1,72 = where a decouples, b loses its kinetic energy and 
becomes a Grassmann Lagrangian multiplier, integrating b out leads to the boundary conditions : 

<i(0) = ( 51 ) 

We also have the trivial boundary conditions 

blM = 6« (0), 6f>(0) = &f f (0) (52) 
Using Eq.|l3|, the above boundary conditions of H' correspond to the following boundary conditions of H: 

«i(0) = -a£(0), 61(0) =6^(0), b s /(Q) = 6^(0) (53) 
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The above boundary conditions can be expressed in terms of the new chiral bosons in Eq.[l9|: 

^lW—^W+tt (54) 

In terms of new physical fcrmions, it reads: 

^i±,i(0) = e ±i *S i±jJl (0) (55) 

This is a NFL fixed point with g — \/2 and the symmetry 0(1) x 0(7) which is the fixed point symmetry of the 
2IK(a) (Fig. fil). At this fixed point, the original electrons scatter into the collective excitations which fit into the S 



spinor representation of SO (8). The finite size spectrum is listed in Table [V 
The local correlation functions at this fixed point are 

(6(r)6(0)) ~ ~ <a s ,,(0, r)a Sii (0, 0)) ~ 4 ( 56 ) 

We can also read the scaling dimensions of the various fields [a] = 0, [b] — 1/2, \b s ,i\ = [b s f] = 1/2, [a s j] = 3/2. 
This NFL fixed point is more unlikely to be observed by experiments, because it has two relevant directions Ai, A_. 
Similar to the discussions in the last section, at the 2IK fixed point, a decouples and b turns into the non-interacting 
scaling field at the fixed point 



The corresponding scaling field in H is 



The impurity spin turns into 



6(T)~ as ,,(0,r) (57) 
-6«(0,r) (58) 



S x (r) ~ iaa s H 

S 9 (r) ~ i(ab s + a s b s ) H 

S z (t) ~ tab s + ■ ■ ■ (59) 
The impurity spin-spin correlation function (S a (T)S a (0)) ~ 1/r. 

On the line of NFL fixed points discussed in the previous section, if 8 = 0, the fixed point symmetry is enlarged 
to the 2IK(b) ( Fig. |l| ). Although these two fixed points have the same symmetry, therefore the same finite size 
spectrum, but the allowed boundary operators are very different. This additional NFL fixed point is also unstable. 
Ai, A_ are two relevant terms with scaling dimension 1/2, any linear combination of the two terms will drive the 
system to a given point on the line of FL fixed points to be discussed in the followiiig section. They will generate two 
dimension 3/2 leading irrelevant operators L_i& s /(0, r), L_i6 s ,i(0, r) respectivelycll and two dimension 2 subleading 

irrelevant operators (energy-momentum tensors ) 6 s /(0, r) dbsl g^ ,T ' 1 , b s i (0, r) abs g^°' r \ As explained in the last section, 
because the two leading irrelevant operators are first order Virasoro descendants, they make no contributions. 72 term 
also has dimension 2. In all there are three dimension 2 subleading irrelevant operators which contribute Ci mp ~ T. 
In all, we have 

C imp ~ T (60) 

From Eqs.|l^,^, we get the susceptibility 

xL P ~9 2 logT (61) 



By Eq.12, the leading and subleading irrelevant operators in H' become 

db 3f (0,T) dV s (0,r) 
dr ' dr 



M°> T ) q- t , 0.(0, r) g^- 

<>, ^db(r) . da'(0, t) 

720(r)-^ - 72 a' s (0, r) ^ ' (62) 
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By EqJT^, the corresponding operators in H are 

db af (0,T) da s (0,T) 
dr ' dr 

6 s /(0,r) , a s (0,r , 7 2 & s (0,r) — 63) 

OT OT OT 



i 



They can be written in terms of the new bosonsl 

— sin$"(0,r), — cos$»(0,r) 
±cos2$:(0,r)-i:(5$';(0,r)) 2 : 

72 (cos2$' s V(0,r) - \ : (0*^(0, r)f :) (64) 

The first order corrections to the single particle L-R Green function due to the 1st and the 2nd operators in Eq.|64| 
vanish. ._. 
Those due to the 3rd and 4th operators areEZI 

(lp 1+ (x 1 ,T 1 )lp\ + (x2,T 2 )) ~ 

J dr(e-i*" (Kl ' ri) (±cos2$^(0,T)- i : (d$ ? s l (0, r)) 2 : ) e i(-*?(-2,r 2 )+^ =Q (65) 
The 5th operator in Eq.^4] makes no contributions to the Green function either: 

( Vi+(zi,ti)VI + (z 2 ,t 2 )) ~ 72 / d r ( e i(-*?^^)+-) e -i*?(^^)) 



x ( e -**"/( :E i. 7 -i)(cos2$ , s 1 / (0,r) - i : (d^ f (0, t)) 2 : ) e i*»/( a:2 '' r3 )) = (66) 
It is easy to see higher order corrections due to the above operators also vanish. 

The y and q terms are two irrelevant operators with scaling dimension 5/2, they can be written in H' as 

: 6(t)^M : M o, r ), : 6(r)^M . M o, T ) (67) 



The bosonized forms in H are 

1 

2 



(cos2$J / (0,r) - i : (9^(0, r)f :)sin^(0,r) 

(cos2^ / (0,r)-i : (9<f^(0,r)) 2 :)cos^(0,r) (68) 



The first order correction due to the first operator isE3 

( ^i+(a:i,ri)Vl + (a;2,7a)>~i / d r ( e -|^(^.n) e i*?Co,r) e 4(-* ? (x 2 ,r 2 )+ 7 r) ) 

x ( e -i*?/(*i.«) : (5*^(0, t)) 3 : e**"< aia '* a >) - - z 2 )-^ 2 (69) 

The above integral is essentially the same as the first part of that in Eq.|35|. 

As explained in the preceding section, due to the factor i difference than the first operator, the 2nd operator makes 
no contribution to the conductivity . ._. 

The other two dimension 5/2 operators ard_3 

L_A/(0,r) = |j^siii$?(0,T), £-20.(0, r) = ^ cos $?(0,r) (70) 

Because they can still be written as total derivatives, therefore make no contributions to the Green functions. 
Conformal transformation of Eq.^ to finite temperature and scaling analysis^ lead to 

a(T) ~ 2<r tt (l + T 3 / 2 ) (71) 

There is no chance to form universal ratios from Eqs j60| , |6l| , [7l| . 

In Appendix C, similar calculations in the old boson basis Eq.|s| are performed. 
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V. THE LINE OF FL FIXED POINTS 



If q — y — 0, this fixed point is located at 71 = 72 = where both a and b lose their kinetic energies and become 
two Grassmann Lagrangian multipliers, integrating them out leads to the following boundary conditions : 

a£i(0) - -a&(0), 6f /ti (0) = -6f /ti (0) (72) 
We also have the trivial boundary conditions 

6^(0)= 6* (0) (73) 

Following the same procedures as those in the discussion of the line of NFL fixed points , the above boundary 
conditions correspond to the boundary conditions in H 

61(0) =%(0), r L (0)=e i9 r R (0) (74) 

The above boundary condition can be expressed in terms of the new bosons in Eq.|l9| 

= $Ir + (75) 

As discussed in the line of NFL fixed points , the physical fcrmions transform as the spinor representation of the 
above boundary conditions, the corresponding boundary conditions are derived in Appendix B, the result is 

= e*"*1& (76) 

It can be checked explicitly the above boundary conditions satisfy all the symmetry requirement. The fermion fields 
with the even and the odd parity under Z2 suffer opposite continuously changing phase shifts ±| along this line of 
fixed points . Depending on the sign of Ao, the impurity will occupy either the even parity state with S x = ^ or 
the odd parity state S x = — i. This simple physical picture should be expected from the starting Hamiltonian EqJ|: 
If we keep Ao term fixed, then V3 term is irrelevant, the y term ( y — 2V2 ) is exact marginal, it causes continuous 
opposite phase shifts for ipi±. V3 term will generate the dimension 2 operator : J z (0)J z (0) :, the OPE of the y term 
and the V3 term will generate the S z J v (0) term, this term will generate another dimension 2 operator : J y (0) J y (0) :. 
The impurity spin-spin correlation function S z (t)S z (0) ~ (7^-) 2 77 which is consistent with Eq.|8l|. 

One particle matrix is Sf — e ±4e / 2 , the residual conductivity <r(0) = 2er u /(l — cos |). This is a line of FL fixed points 
with <7 = 1 and the symmetry U(l) x 0(6) ~ U(A) which interpolates continuously between the non-interacting fixed 
point and the 2CSFK fixed point found by the author in RefEj. If 9 = 0, namely Ai = 0, the residual conductivity 
comes from the potential scattering term discussed in Sec. V which is omitted in this section, the fixed point symmetry 
is enlarged to 0(8) which is the fixed point symmetry of the non-interacting electrons. We can see this very intuitively 
from Eq. |^: if y = Ai = A2 = 0, then they will remain zero, the impurity will be in either the even parity state or the 
odd parity state depending on the sigH of Ao, the electrons and the impurity will asymptotically decouple, therefore 
the electrons are asymptotically freeEJ. Non-zero A 2 will not change the above physical picture, because it will not 
generate y and Ai terms, however if A + = Ao + A2 = 0, the 2IK(b) is reached, if A_ = Ao — A2 = 0, the 2IK(a) 
is reached, As shown in Fig.[l], both fixed points will flow to the free fixed point asymptotically. If 9 = tt, namely 
A_ = 0, <t(0) = 2cr u , the fixed point symmetry is enlarged to 0(2) x 0(6) which is the the fixed point symmetry of 
the 2CSFK. 

This line of FL fixed points is stable. There is no relevant operator. There is one marginal operator in the 
spin sector 9$™(0) which changes the angle 9 in Eq. |7^. 71 and 72 terms which are leading irrelevant operators 
have dimension 2, they lead to typical fermi liquid behaviors, q term has dimension 3; in H' , it can be written as 

: b(r) : 6 S) i(0,r)6 s /(0, r). The y term has dimension 4; in H' , it can be written as : a(r) b(r) :. 
This operator contents are completely consistent with our direct-analysis mentioned above. This complete agreement 
provides a very powerful check on the method developed in RefJ13. 
The local correlation functions at this fixed point are 

(a(r)a(O)) ~ ± (6 s/)i (0,r)6 s/ii (0,0)) ~ % 



T" 
, 2 



<6(r)6(0)) ~ (o a ,,(0, r)a s>i (0, 0)) ~ (77) 

From the above equation, we can also read the scaling dimension of the various fields [a] = [b] = 1/2, [b a< i] = 
1/2, [a a ,i] = [b sfA ] = 3/2. 
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At the fixed point, a,b turn into the non-interacting scaling fields in H' 

o(t) ~ b sfti (Q, t) = J_(-A_6 a)< (0, r) + Aib fl/ (0, r)) 

6(t) ~a s ,*(0,r) (78) 
The corresponding two scaling fields in H are 

--?-(- A_a s + Ai6 s/ ) 
-6.(0, r) (79) 

The impurity spin turns into 

5 x (t) ~ -i-^-a s b s f H 

<Sy( T ) ~ «(a s 6 s + -r— (-A_a s + Ai6 s /)6 S ) H 

Ax 

S»(t) ~i-^-(-A_a s + A!6 s/ )6 s + --- (80) 

The impurity spin-spin correlation function shows typical FL behavior 

(S a {T)S a {0)) ~ 1/t 2 (81) 
The two leading irrelevant operator in H become 

da{r) 



7io(r)- 



<9t 



A 2 da s (r) AAi 8^,,^, gMr) 
1"a — J ° s — 5 l ~\2 — fl s 5 1" I "a — > ° s f a 

= cos6 : cos2$^(0,r) : -- : (<9$ r s l (0, r)) 2 : 
-sin0(: sm2$™(0,r) : -d 2 $"(0,r)) 

72&(r)^ ~ 72O 0082^(0^) : -\ : (d^ f (0,r)f :) (82) 

Although the first operator do depend on the angle 0, its scaling dimension pemains 2, therefore will not affect the 
exponents of any physical measurable quantities. We refer the readers to RefBa for the detailed similar calculations 
on the single particle Green function and the electron conductivity . Second order corrections to the single particle 
Green functions from the two leading irrelevant operators lead to <j(T) ~ <r(0) + C(6)T 2 . 



VI. THE EFFECTS OF EXTERNAL MAGNETIC FIELD 

According to Ref.i, the parameters in Eq.|| 1 



are 



Vi = irppV (83) 
sinfcpi? 

V-2 = ttPfV f 84 
kpR 



w-W 1 -^'" (85> 

The external magnetic field H breaks the SU(2) flavor ( the real spin ) symmetry, it causes the energy band of 
spin | electrons to shift downwards, that of spin j to shift upwards. Channel 1 and 2 electrons have different fermi 
momenta, therefore couple to the impurity with different strength. Setting the external strain h = 0, the Hamiltonian 
is 



14 



H = H +i5v F J dx{^\ a {x) — — tp^ix) — ) 

+ FiJ c (0) + +*Vi Jf (0) + 2y 2 J x (0) + 25^2 J x {0) 
+ W 3 S z J z (0) + A^J^O)^ + J y (0)S v ) 
+ 4SV 3 S Z J Z (0) + 5A 1 (J X (0)S X + J^O)^) 

+ ^is*^ + A 2 27rr c (S-V'I t V ) ilV ) 2T'02J. + (86) 

7TT C 1 

Where J a (x) = Jf (x) — J2 (x) and all the <5 terms are ~ iJ. 

The term 4^- J dx(9$ c (a;)9$/(a;) + 9$ s (x)9 ( i> s /(a;)) does not couple to the impurity, therefore can be neglected. 

It is important to observe the bare hopping term and the two electron assisted hopping term are not affected by 
the magnetic field. Following RefH, the transformed Hamiltonian H' = UHU^ 1 is 

H' = H + 2yaba s>i (0)b sf (0) + 28yaba sf {Q)b s>i {Q) 
+ 2qaba s , i (0)b s ,i(0) + 25qaba sf (0)b s f(0) 

- i-7==ab sf (0) + i _L Sa s/ (0) 

V ZlTT c lltT c 

" A ^=ba s 40) + i-=^=ab Sii (0) (87) 



i i /2tTT c ' \J2"KT C 

Performing the rotation Eq||, the above equation can be rewritten as 

H' = H + 2qa&a S!i (0)b Sti (0) + 2Sqaba sf {0)b sf>i {0) 
+ 2yaba sA (0)b sf 40) + 25yaba sf (0)b Sti (0) 



Ak =a6 s/ii (0) + i-J^6a s>i (0) + i^=ba sf (0) (88) 



\J2-KT C \J2-KT C ' \J2lXT c 

It is evident that the magnetic field H introduces, another relevant operator with scaling dimension 1/2. Under 
the combination of the two relevant directions in the above equation, the system flows to the line of FL fixed points 
with the boundary conditions 

-!>':./. = <Kn + «- , 'Kv./. = >K:j ; + e.f m 

If H = 0, 8 s f = 0, the boundary condition Eq.|75| is recovered. If A + = 0, 6 s f = tt. 

There are two marginal operators along this line of FL fixed points, one <9$"(0) is in the spin sector which changes 
the angle 9 S another 9$™j(0) is in the spin-flavor sector which changes the angle 8 s f. The corresponding boundary 
conditions in the original fcrmion basis can be similarly worked out as in the last section. 



VII. SCALING ANALYSIS OF THE PHYSICAL MEASURABLE QUANTITIES 



In this section, following the methods developed in Ref.caE3 and also considering the correction due to the leading 
irrelevant operators, we derive the scaling functions of the conductivity, impurity specific heat and susceptibility : 

A(T, A + ,H, A) = F(^± %XVT) (90) 

VI VI 

Where a, b are non-universal metric factors which depend on 0, 9 p h and the cutoff of the low energy theory, the 
dependence on 6 is due to the existence of the exactly marginal operator 9$™(0) in the spin sectored, the Kondo 
temperature is given by Tk ~ A~ 2 . 

We confine T < Tk, so \y/T is a small parameter, we expand the right hand side of Eq.^o] in terms of the leading 
irrelevant operator 

A(T, A + ,H, A) = / (^t + AVT/ l( ^, ™) + (XVTf , *§) + • • ■ (91) 
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For simplicity, we only consider A + ^ or H ^ 0. The general case can be discussed along the similar line 

of Ref.EI. 

From Eq.^8|, it is easy to observe that A + term and the magnetic field H term play very similar roles. In the 
following, we only explicitly derive the scaling function in terms of A + . The scaling functions in the presence of H 
can be obtained by replacing A + by H. 

As discussed in Sec. V, depending on the sign of A + , the impurity is either in even parity or odd parity states, but 
the physical measurable quantities should not depend on if the system flows to FL1 (even parity) or FL2 (odd parity 
) , so the above scaling function should only depend on | A + 1 . 

In the following, we discuss the scaling functions of a(T) — cr(0), Ci mp , Ximpi Ximp respectively. Here Ximpi Ximp are 
the impurity hopping and spin susceptibility respectively. 

The scaling function of the conductivity 

We look at the two different limits of the / functions. 

Keeping T < Tk fixed, let A + — > 0, the system is in the Quantum Critical(QC) regime controlled by the line of 
NFL fixed points. We can do perturbative expansions in terms of A + /-\/T. As discussed in Refs.tStS, the overall sign 
ambiguity in the spinor representation of the NFL boundary condition Eq.^0] should be fixed by the requirement of 
symmetry and analyticitya. The perturbative expansions are 

a i+ (T, A + , A = 0) - a(0) ~ % + (%) 3 + • ■ ■ 



T WT 

a,_(T, A + , A = 0) - a(0) ~ - (% 3 + • • • (92) 



The total conductivity tr(T, A + ,A = 0) — ct(0) = 0, therefore fo(x) = 0. The conductivities from the different 
parities have-to cancel each other, otherwise, we get a non-analytic dependence at small magnetic field at finite 
temperatureEJ. 

fi(x) = l + x 2 +x A + ---, x«l (93) 
Substituting the above equation into Eq.|9l|, we get 

AA 2 

a(T,A + ,\)-o-{0) = \VT+^± + --- (94) 

VT 



Keep A_|_ fixed, but small, let T — > 0, the system is in FL1 (or FL2) regime, the conductivity should reduce to the 
FL form 

fo(x) = 

/i(a;) = |x| + |a;|- 3 + ---, x > 1 

Ms) = (95) 

Substituting the above equation into Eq.|9l|, we have 

a(T, A + , A) - a(0) = A|A+| + (A|A+|) 3 + +(A|A+|- 3 + • • -)T 2 + ■■■ (96) 

The above equation indicates that the coefficient of T 2 diverges as A|A + | -3 instead of A^ 4 as we approach to the 
line of NFL fixed points. This means the relevant operator A + with scaling dimension 1/2 combined with the leading 
irrelevant operator A with scaling dimension —1/2 near the line of NFL fixed points will turn into one of the irrelevant 
operators \fl-2 with scaling dimension —2 near the line of FL fixed points 

Afl,-2 ~ A|A + |~ 3 (97) 

First order perturbation in this operator leads to Eq.|96]. 
The scaling function of the impurity specific heat 

In the QC regime, the perturbative expansions give (up to possible logarithm): 

g Q (x) = x 2 + x 4 H , x < 1 

91 (x) = 

92{x) = 1 + x 2 + ■ ••, x<l (98) 
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Substituting the above equation into Eq.|9l], we get 

C imp = ^ + ^ + X 2 TlogT + X 2 A 2 + + --- (99) 

It was known that there are accidental logarithmic violations of scaling whcn-thc number of channel is two!. This 
violation has. nothing to do with the existence of marginally irrelevant operatorscfl. Similar violation occur in itinerant 
magnetismEl 

In the FL regime, the impurity specific heat should reduce to the FL form 

go(x) = x~ 2 + • ■ • , x ~> 1 

92{x) =cH , x>l (100) 

Substituting the above equation into Eq.|9l|, we get 

C imp = T(A- 2 + A 2 + •■■) + •• ■ (101) 

The above equation indicates that the coefficient of T diverges as A^ 2 as we approach to the line of NFL fixed 
points. This means the relevant operator A + with scaling dimension 1/2 near the line of NFL fixed points will turn 
into one of the leading irrelevant operators \fl,-i with scaling dimension —1 near the line of FL fixed points 

Afl,— l ~ a; 2 (102) 



First order perturbation in this operator leads to Eq.101. However, as shown in Eq.|96|, this leading irrelevant 
operator make no contribution to the total conductivity, even though it makes contribution to even and odd parity 
conductivity separately (see Eq.^2| ) . 

The scaling function of the impurity hopping susceptibility 

In the QC regime, the perturbative expansions give (up to possible logarithm) 

X L P = A 2 logT+^± + ... (103) 

In the FL regime 

xL P = A 2 logl/A 2 + + .-. (104) 

The exact crossover function can be calculated along the EK line in Eq.|l0|. In the FL regime, the Wilson Ratio 
R = TXimp/Ctmp ~ \ 2 A 2 + log 1/A 2 + is very small as A+ -> 0. 
The scaling function of the impurity spin susceptibility 
In this part, we set A + = and consider finite H. 

In the QC regime, the perturbative expansions give (up to possible logarithm) 

h (x) = log aTH , x<l 

hi(x) = 

h 2 (x) = c 2 +x 2 H , x<l (105) 

Substituting the above equation into Eq.|9l|, we get 

Xlmp = log aT + \ 2 T + \ 2 H 2 + ■■■ (106) 
In the FL regime, it was shown in RcfE 

Xim P = log aH+--- (107) 
Actually, the whole crossover functions go{x), ho(x) have been calculated along the EK line in RefR 
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VIII. DISCUSSIONS ON EXPERIMENTS AND CONCLUSIONS 




In this paper, we brought about the rich phase diagram of a non- magnetic impurity hopping between two sites 
in a metal (Fig. [l]). As discussed in Sec. IV, the NFL fixed point with the symmetry 2IK(a) is very unlikely to be 
observed, although it has very interesting behaviors Ci mp ~ T, Ximp ~ l°g^ an d °~{T) ~ 2cr u (l — T 3 / 2 ). The peculiar 
behaviors of Ci mp ,Xi m p are due to the 'orbital field' couples to a non conserved current. 

Ralph et al found that the experiment data show y/T behavior for OAK < T < Tki ~ 4K and concluded the two 
sites system fall in the Quantum Critical (QC) regime controlled by the " 2CK fixed point" . They also discussed the 
pver to the FL regime in the presence of magnetic field H which acts as a channel anisotropy of scaling dimension 
and in the presence of asymmetry of the two sites which acts as a local magnetic field of scaling dimension 
As first pointed out by MS, even the two sites are exactly symmetric, therefore, the two channel are exactly 
equivalent, there is another dimension 1/2 operator A + which will drive the system to FL regimen. 

In this paper, we find the " 2CK fixed point " actually is a line of NFL fixed points which interpolates continuously 
between the 2CK and the 2IK(a). As GarO showed that under different canonical transformation than employed in 
this paperEa, the 2IK model can be mapped to the 2CK model. This paper discussed the two apparent different fixed 
points in a unified-Jjt&mework. Although P-H symmetry breaking is a relevant perturbation in the original 2IK model 
discussed in RefsScS, but its effect is trivial in this model. Because the two models have different global symmetry, 
although the fixed point is exactly the same, the allowed boundary operators are different. 

We discovered a marginal operator in the spin sector which is responsible for this line of NFL fixed points. In a real 
system, there is always P-H symmetry breaking, therefore there is always a marginal operator in the charge sector. 
Eq.[l7] show that the coefficient of \/T depend on this breaking any way. However, the marginal operator identified 
in this paper is in the spin sector, it, combined with the leading irrelevant operator which contributes to the \Ff 
behavior of the conductivity, will always generate the dimension 1/2 relevant operator. The existence of the line of 
NFL fixed points and the existence of the relevant operator is closely related. There is no reason why the coefficient 
of this relevant operator is so small that it can be neglected in the scaling analysis in the last section. 

The crossover scale from the weak coupling fixed point q = to any point on the line of NFL fixed points is given 
by the Kondo scale Tki ~ D(Ak ) 2 ~ A~ 2 ( see Eq|x] ), the crossover scale from a given point on the line of NFL 
fixed points to the corresponding point on the line of FL fixed points is set by A ~ D(A + ) 2 , the finite size scaling 
at finite temperature T leads to the universal scaling function Eq. Because there is no reliable way to estimate 
the magnitude of A + which is always present. It is very hard to say if the experimental situation do fall in the QC 
regime controlled by any fixed point on this line. 

The experiment measured the magnetic field dependence of conductance signal. Assuming A + is so small that it 
can be neglected, scaling Eqs |6| in the field H shows that in the FL regime, the conductance should depend on X\H\ 
which is consistent with the experimental findingu. The coefficient of T 2 of the conductivity should scale as H~ 3 . 
Because Tk ~ A ~ 2 ~ 4K", the lowest temperature achieved in the experiment is T m ; n ~ O.IK , if O.IK < \H\ < A -1 ~ 
\/Tk = 2K, Eq.101 shows that the linear T coefficient of the impurity specific heat should scale as H~ 2 . Eq.107 
shows that the impurity susceptibility should scale as log|i?|. So far, there is no experiment test of these scaling 
relations in this range of magnetic field. It should be possible to extract Ximp from experimental data, because the 
impurity does not carriz_real spin. There is no difficulty caused by the conduction electrons and the impurity having 



different Lande factorstfl. 

It is very difficult, but not impossible to measure x\ m p by adding pressure and strain to the system. Because Ximp 
here is the hopping susceptibility. The difficulty caused by the conduction electrons and the impurity having different 
Lande factors is not a issue either. Unfortunately, no universal ratios can be formed among the amplitudes of the 
three quantities except at the 2CK point on this line of NFL fixed points, this is because (1) the strain couples to a 
non conserved current (2) this is a line of NFL fixed points instead of a fixed point. 

The experiment also measured the conductance signal when a finite voltage V is applied to the point contacts and 
find eV/T scaling in the temperature range OAK < T < 4K. It is not clear how the position on this line of NFL fixed 
points enter the expression of the non-equilibrium conductivity calculations. This question deserve to be addressed 
seriousiyEHl. 

Refuel showed that a non-magnetic impurity hopping around 3 sites arranged around a triangular is delocalized and 
the metal shows either the one channel Kondo fermi liquid behavior or the 2CK non-fcrmi liquid behavior. They also 
conjectured that there may be a " NFL fixed point " possessing local conformal symmetry SUs(2) x ST/q^) x U(l) 
separating these two fixed points. The insights gained from the two sites problem discussed in this paper implies 
that "the NFL fixed point" separating the one channel Kondo fixed point and the 2CK fixed point may be a line 
of NFL fixed points instead of a NFL fixed point. The symmetry of this line of NFL fixed points is a interesting 
open question, but it should be smaller than SU(3), just like the symmetry in Xhe spin sector of the line of NFL 
fixed points in the two sites problem is U(V) x 0(1) which is smaller than S't/(2)Eil. The higher symmetry SU(3) is 
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realized at just one point on this line of NFL fixed points . It was shown in RefHa that the 2CK fixed point with the 
symmetry SU(2) can be realized in Cs v or higher symmetry, because it is indeed possible for the ground state of the 
impurity-electrons complex to transform as the La representation of C$ v group, therefore a doublet. This NFL fixed 
point was also shown to be stable. Similarly, ReE3 pointed out that the stable 5f/(3) NFL fixed point can be realized 
in the system of a non-magnetic impurity hopping around the tetrahedral or octahedral sites in a cubic crystal when 
the ground state is a triplet. However, as the symmetry get higher, the NFL fixed point with higher symmetry will 
be more unlikely to be realized by experiments, because the number of relevant processes will increase. 
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TABLE I. The finite size spectrum at the line of NFL fixed points with the symmetry O(l) x U(l) x 0(5) when < 6 < \. 
Eo — jg + §(tt) 2 - NSs is the state achieved by twisting NS by the angle 25, namely = —e l2S ip(l). Rs is the state 

achieved by twisting R by the angle 28, namely ip(—l) = e t2S i/j(l). NS+lst is the first excited state in NS sector et. al. Only 
when J < S < J, the 4th row has lower energy than the 5th row. If 8 = 0, the symmetry is enlarged to O(l) x 0(7), the finite 
size spectrum of the 2IK fixed point is recovered. If 5 = ?, the symmetry is enlarged to 0(3) x 0(5), the finite size spectrum 
of the 2CK fixed point is recovered. 
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TABLE II. The finite size spectrum at the line of NFL fixed points with the symmetry O(l) x [7(1) x 0(5) when — | < 8 < 0. 
Only when — ? < 8 < 0, the 3rd row has lower energy than the 4th row. If 8 = 0, the symmetry is enlarged to O(l) x 0(7), 
the finite size spectrum of the 2IK fixed point is recovered. If 8 = — ^, the symmetry is enlarged to 0(3) x 0(5), the finite size 
spectrum of the 2CK fixed point is recovered. 
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TABLE III. The finite size spectrum of the 2CK fixed point 



19 



o(i) 


0(7) 


-J—(E - —) 


T^porpriprapv" 


R 


NS 





2 


NS 


R 


1 


8 


R 


NS + 1st 


2 


14 


NS+lst 


R 


7 

s 


8 


R 


NS + 2nd 


1 


42 


R + 1st 


NS 


1 


2 



TABLE IV. The finite size spectrum of the 2IK fixed point 
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TABLE V. The finite size spectrum at the line of FL fixed points with the symmetry 17(1) x 0(6) when < S < J. 
i?o = 5 ('I) 2 - K 5 = 0, the symmetry is enlarged to 0(8), the finite size spectrum of the free fermion fixed point is recovered. 
If S = ^, the symmetry is enlarged to 0(2) x 0(6), the finite size spectrum of the 2CSFK is recovered. 
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TABLE VI. The finite size spectrum at the line of FL fixed points with the symmetry U(l) x 0(6) when —5 < 5 < 0. 
So = §(-) 2 ' If £ = 0, the symmetry is enlarged to 0(8), the finite size spectrum of the free fermion fixed point is recovered. 
If 5 = — the symmetry is enlarged to 0(2) x 0(6), the finite size spectrum of the 2CSFK is recovered. 
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TABLE VII. The finite size spectrum of the free fermions with both NS and R sector; 
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TABLE VIII. The finite size spectrum of the 2CSFK fixed poinfl 
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APPENDIX A: THE FINITE SIZE SPECTRUM OF ONE COMPLEX FERMIONS 

In this appendix, we closely follow the notations ofEl. 

The energy momentum tensor of one complex fermion in the complex plane is: 

T(z) =: V* (z)dip(z) : (Al) 
Where the complex fermion can be written in terms of two Majorana fermion: 

iKz) = -j=(M*) + iM*)), 4>*(z) = -j=(M*)-iM*)) (A2) 

The mode expansions of the complex fermion are 

#(z) = ^^z-™- 1 / 2 , V'« = ^=(Vd + #£) 

*vw=X>;*- n - 1/2 > ^ n = ^i-ii>i) (as) 

The modes satisfy the commutation relations: 

{<l>n, VVn} = {iPni V4} = °> {V'n, V4} = <Wm,0 (A4) 

Under the conformal transformation z = e w , mapping the cylinder, labeled by w = t + icr, to the plane, labeled by 

z 

V>H - Yl ^ e ~ nw > V>» = E ^ e_m " (A5) 

n n 

Where n G Z — 6 in order to satisfy the boundary conditions on the cylinder 

ip(w + i2n) = e i27r< V(™), < 6 < 1 (A6) 
The energy of the complex fermion on the cylinder is 

H = L = J2 ( n ~ e ) ■■ ^n-ei>-n+6 ■ (A7) 
nez 

From the commutation relations 

[L , tpn-e] =-(n- 9)ip n -e, [L , ip*_ e ] =-(n- 0)i/j*_ g (A8) 
and the fact that the regularity at r = — oo requires the positive modes annihilate the vacuum 

tp n -e\0)e = 0, n-9>0 (A9) 

We can construct the Hilbert space 

■■■r- m -e^-n-e\0)e (A10) 

Where m, n are no negative integers. 

The energy of the twisted vacuum |0)g is 

Eo^fW^-^ + k-) 2 , e = l+ 6 - (All) 
The energies of the excited states in Eq. Ald| are 

E = E o + m + + n + 0-\ (A12) 

If 8 = 0, then the states are in Ramond ( periodic ) sector, because of zero modes, the ground state degeneracy is 
4, the ground states are |0),V'o|0>,^ol >>V'oV'o|0). If 0=\, the states are in Neveu-Schwarz (anti-periodic) sector. 
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APPENDIX B: THE DERIVATION OF THE BOUNDARY CONDITIONS 



SO (8) has three 8-dimensional representations, one vector representation 8„; two spinor representations, one 8 C 
with positive chirality T = 1, another 8 S with negative chirality r = — 1. They can be bosonized as 
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Where $ c , $ s , are defined by /irsi triality transformation Eq.||. 

^c.s: $ s ,s! $/, s , ^s/,s are defined by second "triality transformations" 

$i T - $ u + $2T + $2j 
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V $i T + $ u - $ 
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Note the only difference between C„ and 5^ is the change of sign of <&ij. 
In the basis of Eq.|5|, we can rewrite Eq.Bl 



( e ~ l ! c \ 



\ 



— i<I> c 



as 



/ 



e -i*2,T 



The other two sets of bosons are expressed in terms of those in EqJ| 



/ *i,T \ 








$ 2;T 




V *2,i / 





$ c - $ s + $/ - $ s/ 





(Bl) 



(B2) 



(B3) 



(B4) 



It is evident that in the first basis of Eq.[5], ip^ transforms as 8 C and C M transforms as 8„, namely, W an d 
exchange roles. It is important to observe that the only difference between tp^ ( T = 1) and S v ( T — — 1 ) is the 
change of sign of $ s . 

Similarly, in the basis of Eq.B2, ip u transforms as 8 S and S v transforms as 8„, namely, ip\ v and S v exchange roles. 
The complex fermions in EqTqare C v fermions. 



1. The derivation by bosonization 



In EqJT|, we bosonize new sets of complex fermions in terms of new sets of bosons $™, $ 
Replacing $ s , $ s / by $™, in Eq.[B4|, we can construct -0™, 5™. 
Eq.|0| and |B4] lead to 

in _ ±10/2 on on _ ±id/2 j n /p>r^ 

where in the exponential, we take +, if a =| and — , if a —I. 

Because the basis [19] is related to the basis || by a 50(8) rotation matrix in 8„, therefore tp" a ,S" a are related to 
"0i,aj5i iQ by the two rotations in 8 C ,8 S respectively. 

We can determine the two rotations in 8 C , 8 S by the following two steps. 

(1) Under $ c — * $ c + A, $/ — * $/ + A, A is an arbitrary angle. 



Eq.B4 dictates 



ipil = c 2 iVi,T + c 22lpt,l (B6) 

(2) Under $ s — > $ s + w, a s — > -a s ,b s — » -& s therefore $™ — > -$™ + tt, — > Under $ s / — > $ s / + n, 

a sf -> -a s/ , 6 s/ -> -& s/ therefore -»■ $^ -> + tt. 
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Eq.B4 dictates 



d = i>i,+ = ^(^*t + V'a) 
= fa,- = "^(^T " V'il) 

Remember S* differ from ip by the change of sign of ^ s , we get the similar relations for S spinor 

S§ — Si,- — —j=[Si-[ — Sii) 



(B7) 



(B8) 



Eqs.B5,B7,B8 give the boundary conditions at the line of NFL fixed points Eq.|21|. 

2. The derivation by 7 matrix 

In this subsection, we follow Ref.0, but use different conventions. 

Jki are the 28 generators of SO (8) and 4 commuting generators W k form a Cartan subalgebra of SO (8) 

(Jkl)mp = SkmSlp ~ ^kp^lm, k,l = 1, • • • , 8 

W k = J 2k -i,2k, fc = l,2,3,4. 



(B9) 



The eight Majorana fermions in Eq.^ transform as 8„. Eq|74| can be written as e ej3S which is a matrix in SO (8) 
(therefore in 8„ ). 

ip u , S u will transform as 16 x 16 spinor representation of 5*0(8). The 16 x 16 matrices can be written in the block 
form 



r = 



7 l 
7*' 



(BIO) 



P satisfy {P,P} = 2S ij if 



rfrfi + yjji' = 25 ij 
y'y + y'y = 2S ij 



The generators of the 16 x 16 representation are 



r« = i[P,P] = J 







7 « 
7 ? 



(Bll) 



(B12) 



Where Yj — H l J = 57 1 '7"' are the generators of 8 C and 8 S respectively. 

Under $ c -> $ c + 9, EqE transform under e eJl2 , from Eq.pl ip^, S M both transform under e e / 2 ( Jl2+j34+j56+j78 ). 
Under <1> S — > $ s + 0, Eq.|]transform under e ej34 , from Eq.B4, ip^jS^ transform under e e ' / 2l y J ^~ J ^+ J ^- J ^) and 
e -e/2{j 12 -j 3i +j 5e -j rs ) reS p ec tively. Under -> $/ + 0, Eq]| transform under e 9Jr ° 6 , from Eq{B4l ^,,5^ botn 
transform under e 9 / 2 ( Jl2+j34 ~ j56 ~' /78 ). Under <f> s / — > $ s / + 0, Eq.|| transform under e eJrs , from Eq.[B4|, ip^, S^ both 
transform under e 6, / 2 ( Jl2+ ' /34 ~ j56_ ' /78 ). Therefore, omitting subscript c, we conclude 



712 = 
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= 7l2,s 
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J 56 
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= -734,s 
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= 778,s 



(B13) 
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The 8 7^ matrices can be expressed as the direct products of 2 x 2 blocks 



7i = e x 1 x (73, 72 = e x 1 x <j\ 

73 = 1 X (73 X e, 74 = 1 X 1 X 1 

75 = e x e x e, 7 6 =(TiX£xl 
77 = 1 x o\ x e, 7s = (73 x e x 1 

Where e = io-i and <7j are the three Pauli matrices. 

It is easy to see 74 is an identity matrix, the other 7 matrices are r eal a nd anti-symmetric. 



(B14) 



It can be checked that indeed 7^ = ^7*7 



3> 7? 



±7 l V satisfy EqjB13 



If Eq.^| tr ansfo rm under e eJa8 , then ^ and S 1 ^ transform under e 6 * 738 and e^ 738 - 3 respectively. 
Using Eq. B14 , straightforward calculation show that e 9738 give Eq.|76|. 

Under the first equation in Eq.[l8], <I> S — > — $ s , from Eq.B4, it is easy to see Vv an d ^ transform to each other, 
therefore lead to Eq.GU 



APPENDIX C: THE CALCULATIONS IN THE OLD BOSON BASIS AT THE ADDITIONAL NFL 

FIXED POINT 

In contrast to the boundary conditions Eqjl^ of the line of NFL fixed points , The boundary conditions Eq.^3] are 
in one of the four Cartan subalgebra, therefore can be expressed in terms of chiral bosons in Eq.|B|: 

$ s , i (0)=-$ s , i j(0)+7T (CI) 

In terms of physical fcrmions, it reads: 

fe,L(0)=e^5 ia ,fl(0) (C2) 

where j a = ± i . 

The five operators in Eq.|6^ can be written in terms of the bosons introduced in Eq.||, 

d d 
— cos$ s/ (0,t), — cos$ s (0,t) 

cos2$ s/ (0,r)-l : (d$ s/ (0, r)) 2 : 

-cos2$ s (0,t)- i : (9$ s (0,r)) 2 : 

7 2 (cos2$ s (0,t) - 1 : (5$ s (0, r)) 2 :) (C3) 

The four dimension 5/2 operators are 

(: cos2$ s (0,t) : -i : (5$ s (0,r)) 2 :) cos $ s/ (0, r) 
(: cos2$ s (0,r) : -I : (d$ s (0,r)) 2 :)cos$ s (0,r) 

J^cos$ s/ (0,t), ^cos$ s (0,r) (C4) 



Very similar arguments to those in Sec. IV lead to <j(T) ~ 2ct„(1 + T 3 ' 2 ). 
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